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Abstract 



We use the relativistic distorted wave impulse approximation to analyze 
data on 16 0(e, e'p) 15 N at \Q 2 \ < 0.4 (GeV/c) 2 that were obtained by differ- 
ent groups and seemed controversial. Results for differential cross-sections, 
response functions and Atl asymmetry are discussed and compared to dif- 
ferent sets of experimental data for proton knockout from p\ m and P3/2 shells 
in 16 O. We compare with a nonrelativistic approach to better identify rel- 
ativistic effects. The present relativistic approach is found to accommodate 
most of the discrepancy between data from different groups, smoothing a long 
standing controversy. 
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I. INTRODUCTION 



Quasielastic (e, e'p) processes are a powerful tool to study bound nucleon properties. 
Indeed, coincidence (e, e'p) measurements at quasielastic kinematics have provided over the 
years detailed information on the energies, momentum distributions and spectroscopic fac- 
tors of bound nucleons. This is so because at quasielastic kinematics the (e, e'p) reaction 
can be treated with confidence in the impulse approximation, i.e., assuming that the de- 
tected knockout proton absorbs the whole momentum (q) and energy (u) of the exchanged 
photon (for recent reviews of the subject see ref. [p]] and references therein). Until recently 
most data were concentrated in the low missing momentum range p m < 300 MeV/c, where 
p m is the recoil momentum of the residual nucleus. In the last years higher p m -regions 
are being probed at small missing energies E m to study further aspects of bound nucleon 
dynamics and nucleon currents. A substantial amount of theoretical work on (e, e'p) has 
been carried out on the basis of nonrelativistic approximations to the nucleon current. This 
is the case of the standard distorted wave impulse approximation (DWIA) ]I[ that uses a 
nonrelativistic approximation to the nucleon current operator and wave functions. DWIA 
has been successfully used over the years [§] to analyze (e, e'p) data using bound and scat- 
tered proton wave functions deduced from phenomenological nonrelativistic potentials. The 
limits of validity of the nonrelativistic DWIA approach are now being studied by Meucci, 
Giusti and Pacati among others. 

In past years we investigated |)|-|| nuclear responses and differential cross-sections for 
exclusive quasielastic electron scattering within the framework of relativistic mean field 
approximations. In the relativistic distorted wave impulse approximation (RDWIA) |5|,p|-[TT[ 
the one-body nucleon current, 

Jjtr(w, = J dpip F (p + g) J^(o;, q)ip B {p) , (1) 

is calculated with relativistic ipB and ipF wave functions for initial bound and final outgoing 
nucleons, respectively, and with relativistic nucleon current operator, J fa. The bound state 
wave function is a four-spinor with well-defined parity and angular momentum quantum 
numbers, and is obtained by solving the Dirac equation with scalar-vector (S-V) potentials 
determined through a Hartree procedure from a relativistic Lagrangian with scalar and 
vector meson terms fI2f . The wave function for the outgoing proton is a solution of the 
Dirac equation containing S-V global optical potentials |13| for a nucleon scattered with 



asymptotic momentum p F . Dirac equations for both scattered and bound wave functions 
are solved in coordinate space and their solutions are then transformed to momentum space 
where necessary. 

Eq. (1) sets up the scenario where differences between RDWIA and DWIA are at play. To 
go from the relativistic to the nonrelativistic approach the one-body (4x4 matrix) current 
operator is first of all expanded in a basis of free nucleon plane waves. This amounts to a 
truncation of the nucleon propagator that ignores negative energy solutions of the free Dirac 
equation. Next, a Pauli reduction is made to transform the current operator into a 2 x 2 
matrix, and an expansion in powers of (q/M) and/or (p/M) (where M is the nucleon mass) 
is made |T4|]. Finally the transition nucleon current is calculated as the matrix element 
between bispinorial, nonrelativistic bound {(f) b) an d scattered {(f) f) wave functions instead 
of the 4-component ipB, ?Pf wave functions. We then cast relativistic effects into: 
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i) Kinematical. These are effects due to the truncation of the current operator to first, 
or higher order in p/M, q/M. For moderate p/M values the relativized form proposed in 

gives proper account of such effects. 

ii) Dynamical. These are effects due to the differences between relativistic and nonrela- 
tivistic wave functions which depend not only on the 4-spinor versus 2-spinor structure, but 
also on the potentials used in the respective Dirac and Schrodinger equations for the bound 
and scattered nucleon. Salient features of dynamical effects are: a) A dynamical depression 
of the upper component of the scattered nucleon wave function in the nuclear interior, typ- 
ically identified as the effect of the Darwin term coming from the derivative of the optical 
S-V potentials 0. b) A dynamical enhancement of the lower components, mainly that of 
the bound nucleon wave function. 

So far, we applied successfully RDWIA to 208 Pb and 40 Ca at low \Q 2 \ [|,|, and to 16 O at 
high \Q 2 \ |]8|Jl5[|. The effect caused by the nonlocal Darwin term for 40 Ca and 208 Pb cases was 
studied in detail in refs. [§,§J. The Darwin term causes an apparent enhanced absorption 
when comparing the RDWIA differential cross-section to the DWIA one at moderate p m 
values, thus predicting larger spectroscopic factors Pfl[T6|j. For larger missing momentum 
values (pm/ (Mc) > 1/3) the lower components of the relativistic wave functions start to play 
a more important role, enhancing the higher momentum components of the nucleon wave 
functions. In previous work J?| we found that RDWIA calculations, compared to standard 
DWIA, tend to produce lower cross sections at p m < 300 MeV/c and larger cross-sections 
at p m > 300 MeV/c, improving agreement |5|-|7j with experiment. 

The effect of the dynamical enhancement of the lower components was studied in RPWIA 
in refs. |[i7|JT8|] . It was also studied in RDWIA in ref. || at high \Q 2 \. In both cases it was 
found to play a crucial role in the TL responses. Recent data (T^| on 16 O at high \Q 2 \ seem 
to confirm former RDWIA predictions. In particular, the richness shown by the structure 
of the Atl asymmetry, which is different for p\/2 and P3/2 shells, is only consistent with 
predictions of relativistic calculations that include the dynamical enhancement of the lower 
components of bound Dirac spinors. Moreover, recent data on polarization observables in 
12 C at \Q 2 \ » 0.5 (GeV/c) 2 also agree nicely with RDWIA analysis US 

For 16 O there is an important controversy in the comparison of theory to data at low 
\Q 2 \- We refer to the data sets from lp-shell proton knockout experiments on 16 performed 
at Saclay [ 2T[] and NIKHEF p2"|j23|| in various kinematics in late 80's/early 90's. These ex- 
periments measured the cross-section as a function of missing momentum and, in particular 
Chinitz et al. J2TJ and Spaltro et al. [ ^2| , extracted also the TL response and Atl asymmetry 
at \Q 2 \ = 0.3 (GeV/c) 2 and 0.2 (GeV/c) 2 , respectively. The measurements from Chinitz et 
al. were compared to relativistic PTJ and nonrelativistic P2j DWIA calculations showing 



relatively small deviations from theory. On the other hand, the data of Spaltro et al. |22 
were compared to results from standard nonrelativistic DWIA calculations, and were found 
to be far from theory. Using nonrelativistic optical potential parameters by Schwandt et 
al. PH , and spectroscopic factors fitted to data in parallel kinematics, Spaltro et al. |22| 
found that the experimental R TL is enhanced by a factor ~ 2.05 for the lp^/2 shell and by 
a factor ~ 1.5 for the lpi/2- 

Though the large discrepancy between DWIA results and experiment found by Spaltro et 
al. may in part be due to two-body currents, calculations of exchange current effects are still 
contradictory |2"5j , |2T)| ] . Hence, the controversy surrounding the TL response and asymmetry 
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data still persists. In view of forthcoming information on 16 responses from experiments at 
Jefferson Lab in the near future, it is important to reexamine these sets of data with RDWIA 
calculations. We investigate whether a systematic fully relativistic analysis of the (e, e'p) 
data at low \Q 2 \ may explain the apparent discrepancies between data from Saclay |H] and 
NIKHEF [0,0]. 

The paper is organized as follows: in Section II we summarize the basic formalism 
needed to describe coincidence electron scattering reactions, paying special attention to 
the relativistic distorted wave impulse approximation (RDWIA). Section III contains the 
theoretical results obtained and their comparison with the experimental data. In section IV 
we present our conclusions. 



II. DESCRIPTION OF (E,E'P) CALCULATIONS 

The general formalism for exclusive electron scattering reactions has been presented in 
detail in several previous papers. We refer in particular to refs. Hl|J§|j27|. Here we just 
summarize the kinematics and focus on those aspects that are of relevance to the points 
under discussion in this paper. As a guide to the reader we write down the unpolarized 
cross-section in Born approximation assuming plane waves for the incoming and outgoing 
electron (treated in the extreme relativistic limit), 

Ka Mott f rec \v L R L + v T R T + v TL R TL cos<p F + v TT R T7 'cos20 F l , (2) 



dVt e de'dVtp 



where e' and Q e are the energy and solid angle corresponding to the scattered electron 
and Qp — {9f,4>f) is the solid angle for the outgoing proton. The factor K is given by 
K = \pf\Ef/ (27t) 3 , with pp the momentum carried by the ejected proton and Ep its energy. 
The term f rec is the usual recoil factor f~J, — |1 — (Ef/Ea-i)(j>a-i - Pf)/\pf\ 2 \, where 
PA-i and Ea-i are the momentum and energy of the residual nucleus, respectively. The 
kinematical factors are vp = A 2 , vp = A/2 + tan 2 # e /2, vpr = A/2, vtl = A y/ A + tan 2 9 e /2 
with A = 1 — (uj/\q\) 2 , where us and q are the energy and momentum transfer in the reaction 
and 8 e the electron scattering angle. The above factors, that contain the dependence on 
the electron kinematics, coincide with those given in [2?|,2"8|j except for a factor v^2 in the 



interference TL term. We remark that in refs. fEILK31 a different convention for K was used 



(see for instance eq. (1) of ref. p^j), which amounts to a factor M/Ep of the responses 
presented in this work with respect to the ones displayed in 0,|2|. 

Our calculation of differential cross-sections and responses includes also the effect of 
Coulomb distortion of the incoming and outgoing electron waves. This breaks the simplicity 
of eq. (0), which is however still useful as a guide. Nevertheless, for 16 Coulomb distortion 
effects in the electron wave functions are tiny (less than 1.5% effect on the cross-section). 

The hadronic current enters only in the response functions R a , a = L,T,TL,TT where 
L and T denote the longitudinal and transverse projections of the nuclear current with 
respect to the momentum transfer q, respectively. Note that the response functions can be 
separated by performing measurements with different kinematical factors v a and/or values 
of the azimuthal angle <pF, while keeping the momentum and energy transfer constant. The 
response R TL is obtained from differential cross-sections at 4>f = 0° and 180°, both in theory 
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and in experiment. Experimental data for the cross-section are often presented in terms of 
reduced cross-sections or effective momentum distributions p{p m ), obtained by integrating 
over a particular missing energy peak the differential cross-section divided by K(2ir) 3 a ep . 
Thus p(p m ) is defined by, 



The free electron-proton cross- sect ion, a ep , is usually taken as occ\ of de Forest PU 



One must be aware that the cross-section given in eq. (g) has a strong dependence in the 
kinematical variables via K and a ep which is removed in the reduced cross-section. For 
instance, at the kinematics of the experiment of Chinitz et al. ^TN (Tp = 160 MeV, \Q 2 \ = 
0.3 (MeV/c) 2 and Ebeam = 580 MeV) a small variation of 5 MeV in T F and u (keeping 
E m and p m constant), may change the cross-section by as much as 7% and the reduced 
cross-section by less than 2%. In order to minimize kinematical dependences, it is safer 
to determine spectroscopic factors by scaling the theory to data on reduced cross sections 
rather than to data on cross-sections. This is so because experimentally, a folding and 
average of the cross-sections, responses and/or reduced cross-sections is performed over 
the experimental acceptance, and central values for the kinematical variables are quoted. 
Theoretical calculations are done for the quoted central values. Due to this, it is not unusual 
that spectroscopic factors may depend on whether one chooses to set the scale by comparing 
to reduced cross-sections or to differential cross-sections, or even to separate responses. In 
this work we first derive the spectroscopic factor (S a ) from the reduced cross-section data. 
Then we use this same factor to compare to data for the individual responses. In this way 
the analysis of R TL and other responses is more consistent and meaningful. 

Another quantity also obtained by the experimentalists and discussed in next section is 
the asymmetry Atl given by 

= q(0 F = 18O°)-q(0 F = O°) 

a(<j) F = 180°) + a{<p F = 0°) ' V ; 

One can see from eq. (2) that this observable is closely related to R TL , with the advantage 
that it is free from the scale factor ambiguity. 



A. Relativistic Distorted Wave Impulse Approximation (RDWIA) 

In RDWIA the process is described H in terms of the one-body nucleon current given in 
eq. ([!]). The relativistic bound nucleon wave function, ipB, is a four-spinor with well defined 
angular momentum quantum numbers K, p, corresponding to the shell under consideration. 
In coordinate space it is given by: 

which is eigenstate of total angular momentum with eigenvalue j = \k\ — 1/2, 

€(r) = J2(^\jP)YUf)xl (6) 

m,cr * 
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with I = k if k > and I = — k — 1 if k < 0. The functions f K , g K satisfy the usual coupled 
linear differential equations [p,30,[3"l"|. 

The wave function for the outgoing proton, ip F , is a scattering solution of the Dirac 
equation, which includes S-V global optical potentials. This wave function is obtained as a 
partial wave expansion in configuration space 0,0: 



ip F (r) = 4tt 1 



IE F + M 
2E F V 



E ■ 



- i5 k l (l m \ a F \j tiYMP^if) 



(?) 



The phase-shifts and radial 



where ip£(r) are four-spinors of the same form as that in eq. ( 
functions are complex because of the complex potential. 

The choice of the current operator J 11 is to some extent arbitrary (see discussion in 
refs. P JT7|j3^| ). Here we consider the two most popular choices denoted as CC1 and CC2 [|29 



fan = Fil" + i^ V Qv 



P + Pf) 



(8) 
(9) 



where F± and F 2 are the nucleon form factors related in the usual way |33| to the electric 
and magnetic Sachs form factors of the dipole form. The variable P in eq. (0) is the 
four-momentum of the initial nucleon for on-shell kinematics, i.e., P M = (E(p),p) (E(p) = 
^/fp + M 2 and p = p F ~q)- 

Thus the evaluation of the one-body current matrix element involves the use of 4 x 4- 
operators and 4-spinors with negative energy components. This is at variance with the 



nonrelativistic (DWIA) approximation where a truncated current operator is used [34] and 
matrix elements are evaluated between bispinorial wave functions (4>b,4>f)- Therefore in 
the discussion of results in next sections we shall refer to relativistic kinematical effects 
-that have to do with the differences due to the use of the complete relativistic current 
operator instead of the truncated one- and to relativistic dynamical effects. A way to fully 
incorporate the kinematical relativistic effects was suggested in ref. ||[14||> an d studied in 
detail in ref. for the reaction 2 H(e, e'p). 



B. Remarks on relativistic dynamical effects 

As mentioned in the introduction the dynamical effects come from the differences between 
relativistic and nonrelativistic potentials and wave functions. In ref. we discussed in 
detail effects on reduced cross-sections for 208 Pb in parallel kinematics due to differences 
between the upper components of the four-spinors if)g, ip F which are Dirac solutions with 
S-V potentials and the bispinors <pB, <Pf which are Schrodinger solutions with standard 
(Wood-Saxon type) potentials for bound and scattered nucleons. 

To illustrate the meaning of this effect we recall that the Dirac equation with S-V po- 
tentials 

{E-f - p ■ ^ - M)?p = (10) 

with 
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E = E — V(r) 
M — M — S(r) 
/ 4>u P \ 

. Ipdown I 
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(11) 

(12) 
(13) 



can be written either as a system of coupled linear differential equations for ip up , ipdown, or 
as a second order differential Schrodinger like-equation for ip up containing also a first order 
derivative term (the Darwin term). Furthermore, using the transformation 



ip up {r) = K(r)<f>(r) 



(14) 



the non-local (Darwin) term can be eliminated to obtain a more standard Schrodinger equa- 
tion with second derivatives only 



2M 



-U 



DEB 



{E 2 -M 2 ) ±/ ^ 

(r) = '-6(f) 

K ' 2M VK ' 



(15) 



with Udeb the Dirac equivalent potential || with central and spin-orbit terms 



Udeb 
V c 

V D 

Vso 
A(r) 



Vc + V S off-i 
' \V 2 - 2EV -S 2 + 2MS + V D 



2M 



1 8 2 A 



1 dA 
rA dr 2A dr 2 
1 1 dA 



3 (dAY 
4A 2 Wfr) 



2M rA dr 
E + M 
E + M 



K 2 (r) 



(16) 



The factor K{r) relating the upper component of the Dirac solution (i/j up ) to the solution 
of the equivalent Schrodinger equation (4>(r)) is called the Darwin factor. As it will be shown 
in next section, K(r) produces a depletion of the outgoing wave function in the nuclear 
interior 

Another dynamical relativistic effect is that coming from the non-zero overlap with Dirac 
sea of the Dirac solutions with S-V potentials. The ipB, i>F wave functions have the general 
structure 



(17) 



where ip^ and ^ > are the projections on the positive and negative energy solutions of the 
Dirac equation for free particles: 



*p {+) (p) = J2 u s(p) u s(p)^(p) = a (+ )(p)v(p) 

s 



(18) 
(19) 
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where we use the notation and conventions of Bjorken and Drell | 3~3|| , so that the positive 
and negative energy projectors are 



A ( ±)(p) 



2M 



(20) 



with P M = (E,p) and E = ^p 2 + M 2 . 

The positive and negative energy components of ip can also be written as 



^~\p) 



s 



(21) 
(22) 



Eqs. (|21| . |22| ) make it more transparent what are the new ingredients of the relativis- 
ts calculation. In particular the difference between ipdown and V4«!m * s what we call the 
dynamical enhancement of the down component. Explicit expressions and figures showing 
a(p) and f3(p) for several orbitals can be found in ref. |18[]. Here we just mention that the 
dynamical enhancement of the down component is proportional to the nonzero Dirac sea 
overlap, /3(p), and that though it is small, it was found to play an important role in the 
TL response function in the RPWIA calculations presented in ref. 
calculations at high \Q 2 \ reported in refs. 



HI, and in the RDWIA 



Its role for the present RDWIA calculations 



is discussed in next section. A way to define an effective 2x2 current operator that includes 
these dynamical relativistic effects was introduced in ref. [Sq . 



C. Projected calculation 

The sensitivity of the different scattering observables to the negative energy components 
can be analyzed by constructing properly normalized 4-spinors of the form in eq. (|18|). Then, 
one can compare the results obtained using the fully relativistic amplitude given in eq. ([I]) 
with those obtained when the negative energy components are projected out. This is done 
when the nucleon current is calculated as 

Jf ++) {u, q) = J dp^\p + q)H", Q)^b\p) , (23) 
where ip^ C0f ) is the positive-energy projection of ips (^f), z.e. ; 

</4 +) (p) = A (+) (^ fl (p), ip£\p + q) = A(+)(p + q)MP + Q) ■ (24) 

The dynamical enhancement of the lower components is contained in the current of 
eq. ([!]), but not in eq. (|23|) . It is important to realize that the positive-energy projec- 
tors inserted in eq. (p3|) depend on the integration variable p. One could also neglect 
this p-dependence by using projection operators corresponding to asymptotic values of 
the momenta, i.e., projectors acting on ifjp and ips respectively, with Pp = (Ef,Pf), 
Pg = Pp — Q^ 1 the asymptotic four-momentum of the outgoing and bound nucleon re- 
spectively, with = (iJ,q) and ZJ = E F — J(pf — q) 2 + M 2 . We refer to this approach 
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as asymptotic projection. This latter projection is almost equivalent to "EMA-noSV" pro- 



cedure employed in ref. fl39| , in which the 4-spinors used have upper components identical 
to the upper components of the Dirac equation solutions, but the lower components are 
obtained with an additional approximation, the effective momentum approach (EMA). Al- 
though EMA-noSV approach also neglects the enhancement of the lower components, it is 
not at all equivalent to the exact projection method in eqs. ( ^3| , p4D . The EMA-noSV ap- 
proach computes the nucleon current with four-spinors that have the same structure than 
the ones encountered in the scattering of free nucleons, because it enforces the relationship 
between upper and lower components to be driven by the asymptotic value of the momenta 
at the nucleon vertex. In particular, the Gordon transformation is exact for EMA-noSV 
approach. Therefore, CC1 and CC2 operators would lead to identical results within EMA- 
noSV, provided the same choices for the off-shell values of u, E, Ep, p and pp are made. 
This would be a strong prerequisite to a factorized calculation, though still not a sufficient 
condition. In order to keep the drawings in section III clear enough, we do not show in the 
figures the results obtained within EMA-noSV, but we shall comment how this approach 
compares with the fully relativistic and/or the projected one. 



III. RESULTS AND DISCUSSION 

In this work we consider three data sets for nucleon knockout from pi/% and p%/2 shells in 
16 that correspond to kinematical conditions of three different experiments. We summarize 
them as follows: 

Set (a) corresponds to the experiment of Leuschner et al. at the Medium Energy Ac- 
celerator (MEA) at NIKHEF-K ||23|| . The coincidence reaction 16 0(e, e'p) 15 N was analyzed 
in quasielastic parallel kinematics at three different beam energies: 304, 456 and 521 MeV. 
The total kinetic energy of the outgoing proton was around 90 MeV. The spectral function 
of 16 was measured in the range, < E m < 40 MeV and —180 < p m < 270 MeV/c, where 
E m and p m are the missing energy and missing momentum, respectively. 

Set (6) corresponds to the experiment performed at the Saclay Linear Accelerator by 
Chinitz et al. ||21|| . The kinematical setup was constant \q\ — to kinematics. The electron 
beam energy was e = 580 MeV, the outgoing proton kinetic energy Tp = 160 MeV, and 
the transfer momentum and energy: |g| = 570 MeV/c and u = 170 MeV {\Q 2 \ = 0.3 
(GeV/c) 2 ). The missing energy resolution was 1.3 MeV, which made not possible to resolve 
the (5/2 + ,l/2 + ) doublet at an excitation energy E x = 5.3 MeV in l5 N from the 3/2~ state 
at E x = 6.3 MeV. 

Set (c), also in \q\ — uj constant kinematics, was obtained by Spaltro et al. [^] with 
the two high-resolution magnetic spectrometers at the medium-energy electron accelerator 
MEA of NIKHEF-K. Data were measured at momentum and energy-transfer values centered 
at (u, \q\) = (90 MeV, 460 MeV/c), i.e., close to the center of the quasielastic peak at 
\Q 2 \ ~ 0.2 (GeV/c) 2 . The experiment covered a missing momentum range from 30 to 190 
MeV/c. The missing energy resolution was about 180 keV, which made it possible to resolve 
the (5/2 + , 1/2+) doublet from the 3/2~ state. 

Next we discuss our results for spectroscopic factors, reduced cross-sections and responses 
corresponding to these sets of data and kinematical conditions. 



S 



In subsection III. A we deduce spectroscopic factors from reduced cross sections, that 
are then used in subsection III.B for response functions. Subsection III. A discusses also 
results corresponding to different relativistic S-V potentials. In our previous work on 40 Ca 
and 208 Pb we found that spectroscopic factors were larger than the ones obtained with the 
nonrelativistic analyses and were very stable when different parameterizations of the S-V 
potentials for bound (HS, NLSH) |,0 and scattered (EDAI, EDAD1, ED AD 2, EDAD3) @ 
protons were used. We shall see that the case of 16 that we examine here is different in 
several respects. 



A. Reduced Cross Section and Spectroscopic Factors 

Let us first discuss the comparison of theory and experiment on reduced cross-sections for 
set (a) (Leuschner et al. ||23|| ). We recall that because of parallel kinematics, for this set the 
only response functions that contribute to the cross-sections are R L and R T . Fig. 1 shows 
the reduced cross-section for pi/ 2 and shells. The sign of p m refers to the projection 
of the initial nucleon momentum along the direction of the transfer momentum q. It is 
defined to be positive for \q\ < \pp\ and negative for \q\ > \pf\- Fully relativistic calculations 
using the CC1 and CC2 current operators (RCC1, RCC2) are shown by thin and thick 
lines respectively. Throughout this paper we use the Coulomb gauge. The Landau gauge 
produces similar results. Gauge ambiguities |4(| are rather small for the fully relativistic 



results in these two gauges HTT|,p~8 



Spectroscopic factors for each of the two shells are evaluated by scaling theoretical calcu- 
lations to experimental data. They are listed in Table I for different choices of wave functions 
and current operators. In this table we also quote the statistical error within parenthesis 
and the \ 2 values per degree of freedom. 

Results on the right panel in Fig. 1 correspond to bound state wave function calculated 
using the parameters of the set NLSH Results with the older HS set |l2|,|3l|, as well as 



with the newest NL3 one E2| are similar. For the scattered proton wave function we use the 



energy-dependent A- independent potential derived by Clark et al. [0 for 16 (EDAI-O). 



Two things are striking in these results that are at variance with the situation we met in 
previous works on 40 Ca and 208 Pb |5]J§: 

1. There are clear deviations in the shape of theoretical and experimental effective mo- 
mentum distributions in the right hand side panel of Fig. 1. Actually, the NLSH wave 
functions have smaller (larger) r.m.s. radii in r-space (p-space) than what is shown 
experimentally. 

2. The spectroscopic factors are small, of the same order or even smaller than nonrel- 
ativistic ones when Perey factor is included in the latter. As seen in Table I the 
spectroscopic factors increase when global A-dependent type (EDAD-1,-2) potentials 
are used instead of the ^-independent potential fitted to 16 (EDAI-O). Moreover, 
the x 2 values are large for NLSH bound wave functions independently of the optical 
potential used. 

We have verified that all EDAD-type calculations (EDAD-1,-2,-3) give similar results 
on reduced cross-sections and responses. Compared to EDAI-O they give about 15 — 20% 
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smaller reduced cross-sections with almost identical shapes. Consequently, EDAD-1,-2,-3 
spectroscopic factors are 15 — 20% larger than EDAI-0 ones but the x 2 values are analogous 
(see Table I). Why this is different from the cases we analyzed in refs. j5],|6j can be easily 
understood from Fig. 2. In this figure we compare the relativistic central potentials (S, V) 
and Darwin factors (K) corresponding to EDAI-0 and to EDAD-1,-2 optical potentials for 
208 Pb (right panels) and for 16 (left panels). We can see that in the case of 16 0, EDAD-1,-2 
potentials produce a deeper K(r), i.e., a larger reduction of the scattered wave in the nuclear 
interior than EDAI-0 potential -also, Vc is somewhat more absorptive-, while in the case 
of 208 Pb both are about the same. Consequently, at the energies considered, EDAD-1,-2 
potential lead to larger spectroscopic factors than EDAI-0 in 16 0, while the two potentials 
lead to similar spectroscopic factors in 208 Pb. The same is true for EDAD-3 and other 
versions of the relativistic EDAD potentials. 

To have a more conclusive determination of 16 spectroscopic factors one would need to 
constrain the optical potential choice by means of inelastic {p,p') data, in addition to the 
elastic ones |43| . But this is not available for the small knock-out proton energies (~ 90 



MeV) considered here nor in a fully relativistic framework. 

The large \ 2 values in the left part of Table I for all the optical potentials have to do with 
the fact that the data do not follow the shape of theoretical reduced cross-sections in the right 
panel of Fig. 1. A similar problem has been found for data sets (b) and (c) where the quality of 
NLSH fits is even worse. This, after all, is not surprising because the standard Lagrangians, 
like NLSH, are fitted to bulk properties of a few heavy nuclei, and one may expect that 
the predicted r.m.s. radii of 16 orbitals differ somewhat from experiment. Unfortunately, 
as seen in Table I this produces large uncertainties in spectroscopic factors. To solve this 
problem we may adjust the parameters of the relativistic potentials (or Lagrangian) so as 
to obtain the correct values of the single particle energies and r.m.s. radii for the orbitals 
considered, in an analogous way to what is usually done in nonrelativistic analyses of (e, e'p) 
data. This is what we do next. 

Compared to data sets (b) and (c), data set (a) has many more data points extending 
over a larger p m range. Therefore, this data set can be used much more reliably to determine 
simultaneously spectroscopic factors and r.m.s. r. values. We have then adopted the following 
strategy: First, we use data set (a) to slightly tune the parameters of the NLSH potential 
so as to reproduce the experimental binding energies and r.m.s.r. values of the pi/2 and 
P3/2 orbitals in 16 0, closely resembling the standard nonrelativistic procedure. We denote 
by NLSH-P the new relativistic potentials and wave functions (see Table II). These new 
relativistic wave functions are then used to make predictions for the kinematical conditions 
of data sets (b) and (c). 

The NLSH-P wave functions are obtained by changing the parameters of the NLSH 
Lagrangian so that the radii and depth of the S and V potential wells derived from the 
Lagrangian are modified in the same proportion. The negative energy content of the resulting 
bound state wave function is barely changed by this procedure. The rescaling of the depth 
size and radii of the NLSH-P wells is within 10% of the initial NLSH ones. The improvement 
obtained in the description of the shapes and quality of the fits is clearly visible in Fig. 1 
and Table I. 

The role played by relativistic dynamical effects is also analyzed from the results pre- 
sented in Fig. 1 and Table I. Each curve in Fig. 1 is scaled by the corresponding spectro- 
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scopic factors in Table I. The reduced cross-sections evaluated after projecting the bound 
and scattered proton wave functions over positive energy states (see eq. ([23D ) are shown 
by thin-dashed (PCC1) and thick-dashed (PCC2) lines. Note that the difference between 
PCC1 and PCC2 results is very small because the so-called Gordon ambiguities are reduced 
after projection fT7| , [T8|| . The results obtained using the asymptotic values of the momenta 



in the projection operator as described in section II. B, are almost identical to the PCC2 
results and thus are not shown here. Once the global scale factor is taken into account, all 
the calculations predict a very similar behavior, what indicates that, aside from the Darwin 
term, the effect of relativistic dynamics in the reduced cross-sections is not important in 
parallel kinematics at low values of \Q 2 \. This agrees with a recent work by Giusti and 



collaborators ||]. This observation also agrees with results of some previous works ||qjl7|jll 
where we saw that the dynamical enhancement of the lower component makes an important 
effect in the cross- section mainly at high missing momentum values and/or in the R TL re- 
sponse function (which does not contribute in parallel kinematics), whereas its influence on 
R L and R T is quite modest. 

Comparing the fully relativistic results with NLSH-P wave functions for CC1 and CC2, 
one observes that the differences are at most of the order of ~ 8%. In the case of the 
projected calculation, we note that the spectroscopic factors are slightly larger than those 
corresponding to a fully relativistic calculation. This is due to the enhancement of the lower 
components of the wave functions which is not contained in the projected approximations. 
Their effect is negligible for CC2 operator and is enhanced by the CC1 choice. We recall 
that another dynamical relativistic effect, namely the Darwin term, is contained in all the 
figures and tables shown here. For EDAI-0 optical potential this effect amounts to a 10% 
reduction of the reduced cross-section in 16 for the kinematics discussed in this work. This 
is comparable to the effect of the Perey factor that was included in nonrelativistic DWEEPY 
calculations |22| , |23| while for EDAD-optical potentials this amounts to a 20% reduction. 

One thus expects the spectroscopic factors listed in Table I for EDAI-0 in the projected 
case to be similar to those obtained from fits with standard nonrelativistic DWIA calculations 
including Perey factors. In this last case the extracted factors for various choices of optical 
potentials are 0.60 < S a < 0.65 for pi/ 2 , and 0.49 < S a < 0.60 for p 3 / 2 which are also 
roughly in agreement with those in Table I for NLSH-P and EDAD-type potentials. 

A smaller spectroscopic factor is expected for the p 3 / 2 shell than for the pi/ 2 , because 
the P3/2 strength is known to be fragmented into three states: the state considered here 
at E m = 18.4 MeV, and two weaker peaks at around 22.0 and 22.7 MeV. According to 
|23| , |4"4|1 the two higher lying peaks would contain about 10% of the total p 3 / 2 strength. The 
spectroscopic factors determined from data set (a) indicate that, taking this extra 10% 
contribution into account, there is similar 3/2~ and l/2~ spectroscopic strength. 

In what follows we use the new bound state wave functions (NLSH-P) to make predictions 
for comparison to the other data sets (b) and (c). We stress that we have used high quality 
data to fix the size of the wave function and that because data set (a) corresponds to parallel 
kinematics no experimental information on the R TL response has been employed. 

Let us now focus on the spectroscopic factors obtained from reduced cross-sections in 
data sets (6) and (c). Fig. 3 shows the reduced cross-sections for py 2 and p 3 / 2 shells. Left 
and right panel correspond respectively to data sets (6) and (c). As in Fig. 1, for each curve 
a global scale factor has been fitted to the experimental data. The corresponding scale 
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factors and their statistical errors are listed in Table III. Similarly to what we saw for set 
(a), also for sets (6) and (c) EDAD-type optical potentials give larger spectroscopic factors 
than EDAI-0 (see also Fig. 4). 

The results for the p%/ 2 shell corresponding to the Saclay experiment (left panel) include 
the contribution of the (5/2 + , l/2 + ) doublet. We have verified that the change in the shape 
of the responses or reduced cross-section after inclusion of the doublet is small. The main 
effect of its inclusion is a decrease of the deduced spectroscopic factor for the p 3 / 2 shell 
of the order of 10%. In ref. [fjl|] the contribution of this doublet was subtracted from the 
experimental data with a procedure based on a nonrelativistic formalism. We have chosen 
to use the uncorrected data from ref. ^1 ] , and include the contribution from the doublet in 
our theoretical calculation. The s — d content has been determined through fits to data set 

(a) for this state P3fl . The values of the spectroscopic factors are £1/2+ = 0.034(2) (RCC1), 
0.034(2) (RCC2), S h/2+ = 0.086(5) (RCC1), 0.088(5) (RCC2) (with a x 2 /dof of the order 
of 0.5). 

Let us now discuss the results corresponding to data set (6) in left panel of Fig. 3. As 
shown, the calculations reproduce in general the experimental data for both shells with the 
scale factors listed in Table III. Although the various approximations give similar results, 
we note that the RCC1 (thin solid line) reduced cross-sections for the pi/ 2 shell are less 
symmetrical around p m = 0, a behavior that is not favored by the data. For this data set 

(b) all the calculations, except RCC1 for the pi/ 2 shell, reproduce well the asymmetry of the 
reduced cross-section. We will return to this point when talking about the TL observables 
in next section. Finally, it is important to remark that the spectroscopic factors obtained 
from the data set (b) (Table III) agree, within statistical errors, with those obtained from 
data set (a) taking into account the systematic error of both experiments: around 5.4% for 
data set (a) |2Bj and 6.3% for data set (6) [fP[ . 



Concerning set (c) |22| the data on reduced cross-sections in right panel of Fig. 3 have 
been obtained from the differential cross-sections and detailed kinematics setup in Appen- 
dices A and D of ref. |45[ (the systematic error for data set (c) reduced cross-section is 



6% f22|). For the pi/ 2 shell, the reduced cross-section is well reproduced by both relativistic 



and projected calculations, except in the case of the RCC1 calculation (thin solid line) that 
underestimates the data for negative missing momentum values. This is consistent with the 
results previously discussed for data set (b). For this shell, the spectroscopic factors that 
fit data set (c) are larger than the ones derived from data sets (a) or (b), but they are all 
compatible within statistical errors. In the case of the P3/2 shell, although the shape of the 
cross-section is well reproduced by the various calculations, the situation on the spectro- 
scopic factors is clearly different (see Table III). With EDAI-0 optical potential the values 
of the spectroscopic factors that fit the p^/ 2 data on reduced cross-sections in set (c) are 
25 — 30% larger than the ones obtained from data set (a). These scale factors are also larger 
than the ones obtained from data set (b), but in this case the discrepancy is of the order of 
15%, which is comparable to the combined systematic and statistical error for these values. 
EDAD-type potentials not only give larger spectroscopic factors but also give, on average, 
better agreement between p 3 / 2 spectroscopic factors of the three different sets (a), (6) and 
(c). This is seen in detail in Tables I and III and is further illustrated in Fig. 4. 

In summary, the shapes of the reduced cross-sections are well described by all the RCC2 
calculations and data sets, what makes us conclude that we can rely on the spectroscopic 
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factors derived with EDAD-1 and NLSH-P potentials. Thus the differences in the P3/2 spec- 
troscopic factors (see Fig. 4) obtained with the same ingredients (wave functions, operators 
and optical potentials) may be attributed either to a global scale variation among the three 
experiments for the p%/2 shell, or to limitations of the theory. Coupled channel contributions 
or MEC could possibly make a different effect for the three kinematics analyzed in this work. 



B. Response Functions and longitudinal-transverse Asymmetry 

In this section we present results for the response functions and asymmetries and com- 



pare them to the data in sets (b) and (c) measured at Saclay jT| and NIKHEF p2] , respec- 
tively. As already mentioned, these two experiments were performed under |g| — u constant 
kinematics so that the TL response and asymmetry (R TL , Atl) can be obtained from the 
cross-sections measured at (pp = 0° and (pp = 180° with the other variables (uj, Q 2 , E m , p m ) 

q 2 

held constant. Moreover, the response functions R L + ~^R TT and R T were also determined 



for data set (c) |22 



Figs. 5 and 6 show respectively R TL and A TL for p 1 / 2 (upper panels) and p 3 / 2 (lower 
panels) corresponding to set (b) (left panels) and set (c) (right panels). In each panel 
we present four curves with the same conventions as in previous figures: RCC1 (thin solid), 
RCC2 (thick solid), PCC1 (thin-dashed) and PCC2 (thick-dashed). Each R TL curve is scaled 
with the corresponding spectroscopic factor quoted in Table III. As it was also the case for 
reduced cross-sections, there are no appreciable differences in the shapes of curves obtained 
with the different types of optical potentials. Obviously, the asymmetry Atl is independent 
on the value of the spectroscopic factor. The results for the p%/2 shell in bottom- left panel of 
Figs. 5 and 6 include the contribution of the (5/2 + ,l/2 + ) doublet as explained for set (b) in 
previous section. The asymmetry Atl was not produced by the Saclay experiment (set (6)), 
but we have deduced Atl from the data using the R TL values as well as the cross-section 



data in 

In Fig. 5 we also show by dotted lines the nonrelativistic results of ref. |22| . For compar- 
ison to previous studies in refs. [^1],^2|], we quote in Table IV the factor required to scale the 
theoretical predictions to the R TL data, additional to the factors in Table III. A value of one 
in this table indicates that the same spectroscopic factor fits both the reduced cross-section 
and R TL , i.e., indicates that the TL strength is consistently predicted by the theory. 

Let us first discuss the comparison between theory and experiment for data set (b). 
From the results shown in Figs. 5 and 6, it is clear that the effects of the negative-energy 
components show up more in R TL and Atl than in the cross-sections (Fig. 3). In the case 
of the pi/2 shell (left-top panel of Fig. 5), the RCC2 calculation agrees with experimental 
data within statistical errors, while PCC1 and PCC2 results for R TL (dashed lines) lie 
about a 30-50% below the data, and the RCC1 calculation (thin solid line) overestimates 
the R TL response by around 20% (see Table IV). In the case of the p 3 / 2 orbit (left-bottom 
panel), all the approximations predict similar curves: The projected results are much closer 
to fully relativistic ones than for the pi/2 shell. Overall, the fully relativistic calculations 
seem to be favored by the data. The fact that in this shell the variation introduced by 
the negative energy components is much smaller than for the pi/2 shell explains why the 
difference between RCC1 and RCC2 results is smaller for the P3/2 than for the pi/2 shell. 
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These results agree with the conclusion reached from RPWIA calculations in ref. |T8[ about 
the behavior of j = £ ± 1/2 spin-orbit partners which was also corroborated in RDWIA 
calculations at high \Q 2 \ ||. 

With regards to the TL observable, independent on the spectroscopic factor, we may 
conclude that for pi/ 2 shell, Atl is best reproduced by RCC2 results, while for p 3 / 2 shell the 
four theoretical results are very close together, and the experimental data agree with all of 
them. 

In the right panels of Figs. 5 and 6 we see the results corresponding to data set (c). 
Most of the comments on data set (b) apply also here, though the data are somewhat more 
scattered and have larger error bars. In the case of the pi/% shell, PCC1 and PCC2 results 
are very similar and lie below the data; among the fully relativistic calculation, the RCC2 
result reproduces the data within statistical errors, while RCC1 overestimates them by a 
35%. In the case of the p 3 / 2 shell (bottom-right panel), all the calculations underestimate 
the experimental TL response by around 17-28%, except RCC1 for which the "additional" 
factor in Table IV is compatible with one within statistical errors. 

In Fig. 7 we show the results for the responses R L + vtt/vlR tt (top panels) and R T 
(bottom panels) for the pi/2 and P3/2 shells compared to the data from NIKHEF [p2]| . Each 
curve is scaled with the spectroscopic factors quoted in Table III. Notice that these responses 
are rather insensitive to dynamical enhancement of lower components. This is consistent 
with the behavior observed in Fig. 1 and also with results of RPWIA calculations |1| . The 
results in Fig. 7 indicate that the separated responses are in general well reproduced by the 
relativistic as well as by the projected calculations for both shells, exception made of the 
data point at the lowest missing momenta where, as indicated by the large error bars, the 
L/T separation is more problematic. 

Summarizing, for the pi/2 shell the RCC2 results agree well with all observables and data 
sets, while RCC1 (projected) calculations show a too large (small) R TL and Atl- For the P3/2 
shell the theoretical calculations lie much closer together, and generally agree with all data 
sets and observables, except for R TL and Atl of data set (c). Although the R TL , Atl data 
on p 3 / 2 in set (c) lie higher than theoretical calculations, they are almost compatible with 
RCC1 and RCC2 calculations within statistical errors. This situation is quite different from 
the one found in ref. |22]| , which is also shown for comparison in Fig. 5. The dotted lines in 
this figure show the nonrelativistic results of ref. that were obtained with nonrelativistic 
spectroscopic factors (0.61(3) for py 2 and 0.53(3) forp 3 / 2 ) and standard (Woods-Saxon type) 
nonrelativistic optical potentials and bound wave functions. The latter were also fitted to 
Leuschner data p3| . 



C. Further Discussion on TL observables 

In this section we focus in Fig. 5 comparing our results to previous nonrelativistic ones. 
The dotted curves -representing the nonrelativistic calculations by Spaltro et. al. p2fl - clearly 



underestimate TL responses for all shell and data sets. The deviation from data is larger 
for p 3 /2 shell, particularly in data set (c), where the dotted curve gives roughly one half of 
the experimental TL response. Why is it that relativistic results in this figure are so much 
closer to data that nonrelativistic ones?. 
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We have examined in detail effects due to the various aspects that are relevant in compar- 
ing relativistic to nonrelativistic results. The effects of Darwin term are already taken into 
account as they basically affect the spectroscopic factors. The effects of the negative energy 
components, as already mentioned, are very small for R L , R T responses in all data sets (a), 

(b) and (c) but, as seen in Fig. 5 and Table IV, they are important for TL observables in 
data sets (b) and (c), particularly for the pi/ 2 shell. We are then left to consider the effect 
of truncation of the current operator (TCO). TCO produces also a negligible effect at the 
kinematics of data set (a), but it is more important at the kinematics of data sets (b) and 

(c) . This again affects more the TL responses and asymmetries where it may represent up 
to a 15% effect (see also ref. ||]). Thus for py 2 shell, TCO roughly explains the difference 
at the maxima between dotted curves and the curves obtained with projected calculations. 
However for p 3 / 2 shell TCO explains only a small fraction of the difference between dotted 
curves and results of projected calculations. The largest fraction of this difference is due to 
the use of a too small spectroscopic factor (see ref. that was taken from ref. and 
that by no means fits the data on reduced cross-sections in set (c). As seen in Tables I, III 
and Fig. 4 the spectroscopic factor deduced from reduced cross-sections in data set (c) is 
25 — 30% larger than that from data set (a). 

The message from this is, not only that relativistic effects are important in perpendicular 
kinematics at low \Q 2 \, but also that a careful analysis of all pieces of information has to 
be done to get a consistent picture of the three different sets of data. Since R TL responses 
are known to be sensitive not only to relativistic effects but also to exchange currents, or 
other possible many-body effects, it is important to establish a clear framework that allows 
to look for the proper magnitude of such effects. 

Indeed if we compare our results to data for the A?l observable that is free from spec- 
troscopic factor ambiguities, we find that all data are well reproduced with the standard 
CC2 current operator, except the p 3 / 2 data in set (c) which are only larger than theory by a 
factor ~ 1.17. This is to be compared to the 2.05 factor that one could expect from ref. ||22|| . 

IV. SUMMARY AND CONCLUSIONS 

In summary we find that the fully relativistic treatment improves substantially the de- 
scription of reduced cross-sections and individual responses of all three sets of data on 
16 0(e,e'p) at low \Q 2 \. Although predictions from CC1 and CC2 current operators are 
rather close in most cases, data seem to favor the CC2 current operator. Therefore our 
remarks here will focus mainly on results with CC2 and with the improved NLSH-P bound 
nucleon wave functions, that have the correct r.m.s. radius. Using the most complete set 
of data on reduced cross-sections in parallel kinematics of Leuschner et al. |2l| (set (a)) we 
obtain spectroscopic factors ranging from 0.58 to 0.64 for pi/ 2 and from 0.45 to 0.55 for 
P3/2, depending on whether we use A-independent (EDAI-O) or A-dependent (EDAD-1,2,3) 
optical potentials. In 16 O, the latter potentials produce a larger Darwin effect, thus larger 
spectroscopic factors. Compared to the cases studied in previous works on 40 Ca and 208 Pb, 
the determination of spectroscopic factors in 16 O with the relativistic approach is different 
in several respects. In the former cases, the standard NLSH wave functions were found to 
reproduce well the shapes of reduced cross-sections and the only fitted parameter was the 
spectroscopic factor. The latter was practically independent on the optical potential used 
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and was ~ 0.7 for the levels just below the Fermi level. On theoretical grounds smaller 
spectroscopic factors for 16 O are expected. In particular, from shell model Monte Carlo 
calculations on 16 |[46||, one may expect S a ~ 0.5 though other theories predict somewhat 



larger values ||7||8| . Larger spectroscopic factors are obtained from Spaltro et al. data [22 



on reduced cross-sections in perpendicular kinematics (set (c)), while Chinitz's et al. data 



21| also in perpendicular kinematics (set (b)) give similar spectroscopic factors than set (a) 



As one can see in Fig. 4, within error bars, spectroscopic factors derived from all data sets 
with EDAD-1 are compatible with each other. To overcome the uncertainty due to the op- 
tical potential (see also Fig. 4) one would need to fit the relativistic potential to both elastic 
and inelastic proton scattering data from 16 in a manner similar to what has been done 
for nonrelativistic potentials ||43|| . The analyses of individual responses is practically inde- 
pendent on the optical potential, once they are scaled by the corresponding spectroscopic 
factors. 

There is a long standing controversy surrounding the TL data for the p\/ 2 and P3/2 shells 
measured at Saclay pi and NIKHEF [22|. We have therefore paid particular attention to 



TL responses and asymmetries and we conclude that there is not a fundamental inconsis- 
tency. Even at the low \Q 2 \ values considered here, the TL response is much more sensitive 
than L and T responses to relativistic effects, in particular to the dynamical enhancement 
of the lower components. The role played by the latter is appreciated comparing fully rel- 
ativistic results (RCC2 or RCC1) to those obtained using wave functions projected on the 
positive energy sector (PCC2 or PCC1). RCC2 results agree well with experimental TL 
responses on pi/ 2 (as well as with TL asymmetries) which are underestimated by PCC2 and 
overestimated by RCC1, because CC1 current operator overemphasizes the role of negative 
energy components. The overall agreement with data on TL responses and asymmetries 
from set (6) and set (c) is quite satisfactory, with the exception of data on p 3 / 2 shell from 
set (c), but even in this case theory is much closer to experiment than previously found 
in ref. |]22|| . In particular, the large difference between data on TL responses from the two 
different sets is well accounted for by the present analyses. This is in contrast with the 



situation depicted in ref. [22], which is represented by dotted lines in Fig. 5. 



In short, the puzzle of the large discrepancy in the TL- response obtained in Saclay |2T 
and NIKHEF |22] , and the "additional" TL-strength found in both experiments is, to a large 
extent, explained by the effect of the negative energy components in the wave functions -a 
dynamical relativistic effect that may not have been expected at low transfer and missing 
momentum. 

The large general mismatch of data set (c) on p 3 / 2 shell seems to point to a normalization 
problem which would require experimental verification. Our analyses indicate that the 
problem is not so much connected to the TL response, but rather to the normalization used. 
Nevertheless, since meson exchange currents and particularly isobar currents are claimed to 
affect more the p 3 / 2 than the pi/ 2 orbitals [l"9"|, it would be interesting to see whether our 
fully relativistic calculation extended to include the isobar and other meson exchange effects 
would lead to better agreement with TL p 3 / 2 data from set (c). It will also be interesting to 
see how relativistic and nonrelativistic approaches compare to new data expected from future 
experiments that have been approved to measure reduced cross-sections and TL-responses 
in 16 with unprecedented precision at Jefferson Lab. 
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TABLES 







NLSH 


NLSH-P 






Pl/2 

CC1 CC2 


P3/2 

CC1 CC2 


Pl/2 

CC1 CC2 


P3/2 

CC1 CC2 






RELATIVISTIC 


EDAI-0 


Q 

°at 

X 2 /dof 


0.58(1) 0.64(2) 
6.6 4.5 


0.45(3) 0.49(3) 
25.3 15.7 


0.54(1) 0.58(1) 
1.3 1.3 


0.43(1) 0.45(1) 
2.7 3.5 


F/DAD-1 

1 J 1 J ii 1 J ± 


o a 
X 2 /dof 


0.63(4) 0.72(2) 
9.6 3.7 


0.56(3) 0.62(2) 
15 7.2 


0.58(1) 0.64(1) 
1.2 1.2 


0.52(1) 0.55(2) 
1.2 4.8 


F/DAD-9 


f 

"a 


0.61(4) 0.69(3) 
10 2.6 


0.53(3) 0.59(2) 
18 9.2 


0.56(1) 0.62(1) 
1.4 1.1 


0.50(1) 0.52(1) 
1.7 4.1 






PROJECTED 


EDAI-0 




0.65(2) 0.66(2) 
4.5 3.2 


0.51(3) 0.52(3) 
16.9 13.3 


0.58(1) 0.59(1) 
1.2 1.6 


0.47(1) 0.46(1) 
3.6 4.3 


EDAD-1 


S a 
X 2 /dof 


0.72(3) 0.74(2) 
4.0 2.6 


0.64(3) 0.64(2) 
7.6 5.9 


0.64(1) 0.65 (2) 
1.4 1.6 


0.57(2) 0.56(2) 
4.4 6.3 


EDAD-2 


S a 
X 2 /dof 


0.69(3) 0.71(2) 
4.6 3.1 


0.61(3) 0.61(3) 
9.7 7.4 


0.62(1) 0.63(1) 
1.3 1.5 


0.55(2) 0.54(2) 
3.7 5.3 



TABLE I. Spectroscopic factors derived from Leuscfmer's experimental reduced cross sections 
in ref. [^3] (data set (a)) using NLSH and NLSH-P relativistic bound nucleon wave functions, and 
EDAI-0, EDAD-1 and EDAD-2 relativistic optical potential parameterizations (see text). Results 
with EDAD-3 are almost identical to the ones with EDAD-1. The numbers within parentheses 
indicate the statistical error. 





Pl/2 


P3/2 




b.e. (MeV) 


r.m.s.r.-r (fm) 


r.m.s.r.-p (MeV) 


b.e. (MeV) 


r.m.s.r.-r (fnL 


r.m.s.r.-p (MeV) 


NLSH 


11.4 


2.838 


175.7 


18.8 


2.679 


185.2 


NLSH-P 


12.1 


3.043 


170.6 


18.4 


2.907 


173.6 



TABLE II. Comparison of binding energies and r.m.s.r. radius in p— and r— space for the 



wave functions NLSH [41] and NLSH-P. The contribution from the negative energy components to 
the norm of the wave function is about 2% in all cases. 
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Set (b) Chinitz et al. |21] 


Set (c) Spaltro et al. @ 




Pl/2 

CCl CC2 


P3/2 

CCl CC2 


Pl/2 

CCl CC2 


P3/2 

CCl CC2 


EDAI-0 (R) 
EDAI-0 (P) 


0.54(4) 0.56(3) 
0.59(4) 0.56(4) 


0.49(2) 0.51(2) 
0.53(4) 0.53(3) 


0.57(3) 0.61(2) 
0.66(2) 0.63(2) 


0.56(1) 0.59(2) 
0.61(1) 0.61(2) 


EDAD-1 (R) 
EDAD-1 (P) 


0.59(4) 0.61(3) 
0.65(4) 0.62(3) 


0.53(3) 0.55(3) 
0.57(5) 0.57(4) 


0.68(4) 0.72(2) 
0.79(3) 0.74(3) 


0.62(2) 0.67(2) 
0.69(3) 0.69(3) 



TABLE III. Spectroscopic factors derived from two different sets of data on experimental re- 
duced cross sections from the full relativistic approach (R) and from the projected one (P). The 
nomenclature used is the same as in Table I. The numbers within parentheses show the statistical 
error only. All results correspond to the NLSH-P bound wave function. 





Set (b) 




Set ( 


c) 






Pl/2 




P3/2 




Pl/2 




P3/2 






N TL 


X 2 /dof 


N TL 


X 2 /dof 




X 2 /dof 




X 2 /dof 


RCC1 


0.83(10) 


0.65 


0.95(17) 


5.3 


0.63(10) 


1.1 


1.09(12) 


3.2 


PCC1 


1.32(42) 


4.5 


1.14(22) 


5.8 


1.15(17) 


1.1 


1.28(12) 


2.5 


RCC2 


1.14(12) 


0.49 


1.02(15) 


3.4 


0.90(13) 


1.0 


1.17(12) 


2.5 


PCC2 


1.48(32) 


2.2 


1.11(17) 


3.8 


1.26(19) 


1.0 


1.28(11) 


2.1 


NR 


1.56(12) 




1.66(9) 




1.50(12) 




2.05(10) 





TABLE IV. Extra scale factor Ntl needed to fit the experimental R response. These factors 
would multiply those in Table III to scale theory to experiment on R TL . A value of one indicates 
that no extra enhancement or quenching of the response is found. The numbers in parentheses 
show the statistical error only. The quality of the fit (x 2 /dof) is also quoted in every case. NR 
corresponds to the nonrelativistic analysis of ref. |2^| . These numbers correspond to EDAI-0 
potential. Very similar numbers are obtained with EDAD-1, EDAD-2 or EDAD-3. 
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FIGURES 



NLSH-P NLSH 

100 p 1 1 1 1 1 1 1 1 1 q p 1 1 1 1 1 




Missing momentum (MeV/c) 



FIG. 1. Reduced cross-sections for proton knockout from lp\/2 and lpz/2 orbits in O versus 
missing momentum p m corresponding to the experiment performed by Leuschner et al. |23| (set (a)). 
The bound relativistic proton wave function has been obtained with the NLSH (right panel) and 
NLSH-P (left panel) parameterization. Theoretical results shown correspond to a fully relativistic 
calculation using the Coulomb gauge and current operators RCC1 (thin solid line) and RCC2 (thick 
solid line). The optical potential used is EDAI-0 from ref. |ll|. Also shown are the results after 
projecting the bound and scattered proton wave functions over positive-energy states: PCC1 (thin 
dashed line), PCC2 (thick dashed line). EMA-noSV results (not shown) are practically identical 
to PCC2 ones. Each curve is scaled by the corresponding spectroscopic factor in Table I. 
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FIG. 2. Real and imaginary part of the optical potentials (upper panels) and real part of the 
Darwin factor (lower panels) (the imaginary part is negligible) for 16 (left panels) and 208 Pb 
(right panels). 
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FIG. 3. Same as Fig. 1 for the experiments performed by Chinitz et al. |2lJ (left panel, set 
(6)) and by Spaltro et al. p2[ (right panel, set (c)). In all the cases the NLSH-P relativistic bound 
proton wave function and EDAI-0 optical potential has been used. For p 3 / 2 shell in set (b) the 
contribution from the nearby 5/2 + and l/2 + states has been taken into account (see text). Each 
curve is scaled by the corresponding spectroscopic factor in Table III. 
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FIG. 4. Spectroscopic factors derived within the fully relativistic approach from the low-Q 2 
data discussed in this work with NLSH-P wave function, CC2 current operator and EDAI-0 (left) 
or EDAD-1 (right) optical potentials. The inner error bars include statistical errors only, the outer 
one includes also the additional systematic error in the reduced cross-sections for each experiment. 
The bands covering the whole \Q 2 \ range corresponds to the value obtained from the data set (a) 
[2J§, while the dots at \Q 2 \ = 0.2 (GeV/c) 2 and 0.3 (GeV/c) 2 correspond to the data set (c) 
and set (b) j2lj, respectively. 
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FIG. 5. Response R TL for proton knockout from ie O for lpi/2 (top panels) and lps/2 (bottom 
panels). Results and data shown correspond to kinematics of data set (b) |2l|] (left) and set (c) |22|| 
(right). Line conventions as in Figs. 1 and 3 (NLSH-P wave function and EDAI-0 optical potential). 
The curves have been scaled by the spectroscopic factors in Table III. Additional dotted curves 
correspond to the nonrelativistic analyses of ref. 22]. 
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FIG. 6. Same as Fig. 5 for the Atl asymmetry. We recall that this observable is independent 
on the spectroscopic factor. 
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FIG. 7. Response functions R L + vtt/vlR TT and R T for the kinematics of data set (c) p2fl . 
Curves and calculations as in Fig. 3. The theoretical results are scaled with spectroscopic factors 
for this same experiment in Table III. 
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